Low-energy magnetic defects in Fe-based superconductors 
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In Fe-pnictides, structural and magnetic planar defects proliferate over differing length scales, 
affecting descriptions of magnetism and superconductivity. "Ordered" Fe moments (0.8-1.04 jib) 
assessed by neutron scattering in the antiferromagnetic state are half the ~1.6 \±b from density- 
functional theory (DFT), while local probes are closer to DFT - a puzzle not yet understood. We 
study via DFT the stability and magnetic properties of competing antiphase and domain boundaries, 
twins, and isolated ncmotwins (equivalent to a twin nuclei), which exhibit low moments confined 
near defect boundaries. A single local- moment picture is thus inappropriate for long-range magnetic 
order. The nanoscale defects are very low energy (starting at 22 raJ/m 2 ), but twins remain so at the 
mesoscale, as observed, and our estimated distances between twins coincide with the observed mag- 
netic correlation length. iVanotwins also explain features in measured pair distribution functions. 
All these defects can be weakly mobile and/or have fluctuations that will lower assessed "ordered" 
moments from longer spatial and/or time averaging. 

PACS numbers: 74.20.-z, 74.25.Ha, 75.25.-j, 75.30.Kz 



Fe-based superconductors (FeSCs) are providing new 
means to understand unconventional superconductivity. 
Due to its ease of synthesis, BaFe2As2 has been chosen as 
a prototype for these systems, where its low-temperature 
(T < 140 K) ground state is a striped, antiferromag- 
netic (AFM) orthorhombic (Fmmm) structure, 1 often 
called a spin-density wave (SDW), and which is repro- 
duced in DFT calculations. At Neel T^ (140 K), both a 
magnetic and structural transition occurs to a tetragonal 
(14/mmm) paramagnet. 1 By doping with a transition- 
metal on the Fe-site, superconductivity (SC) can be 
achieved, and similarly with pressure^™ There are strong 
connections between the magnetism and SC. Dopants 
weaken the magnetic state and Cooper pairing is, per- 
haps, driven by increased magnetic fluctuations out of 
the ground state. 5 6 Fermi-surface (FS) nesting is ap- 
parent from DFT calculations and agrees with angle- 
resolved photoemission (ARPES), suggesting an itiner- 
ant nature^S and which is supported from the spin- 
wave dispersion.^ Furthermore, DFT explains quanti- 
tatively effects of doping on FS nesting, and why Cu 
doping behaves differently than Co and NiP^ K £C Fe2Se2 
(isostructural to BaFe2As2) does not have the hole pock- 
ets needed for FS nesting,^ as DFT finds Also, local, 
multi-orbital models yield a correct groundstate with a 
s ± pairing symmetry for the gap function! 14 * 15 ^ So, mag- 
netism may be a result of itinerant or local electronic 
interactions with the likely description in between, in- 
dicated by measured fluctuating moments (m 2 ) « 3.2 fi% 
that are much less than that given by Hund's ruleP^ 
Our analysis of magnetic defects also supports this con- 
clusion, where we find locality and multiple magnetic 
states at defect boundaries. DFT results for BaFe2As2 
also show a strong coupling between the structure and 
magnetism!^^ We explore this by varying the structural 
parameters in the presence of magnetic defects. 



While DFT supports the observed SDW for th e parent 
compound, the Fe moment (1.6 — 1.9 /iB / 19 | 2Q| is twice 
that assessed for the average "ordered" moment from 
neutron diffraction (0.8 - 1.04 fi B )^^ In fact, vari- 
ous experiments assess very different Fe moments. Core- 
electron spectroscopy^ finds 2.1 like DFT, while 
57 Fe Mossbauei^ and nuclear magnetic resonance^ find 

0. 81 he, as in diffraction assessments. For Fe-based mag- 
nets such a large discrepancy between ordered moments 
from theory and experiment is unusual. DFT+DMFT 
results show that there is large range in values of Hub- 
bard U (0.5-4 eV) depending on the number of bands 
included^ nonetheless, due to Fe(3d)-As(4p) hybridiza- 
tion, the electronic structure is weakly correlated 16 2 ^and 
this agrees with x-ray and photoemission spectra. 26 Spin- 
orbit and hybridization (controlled by Fe/As planar spac- 
ing) in a U-based model was used to explain the small 
in-plane moments in Fe-pnictides Yet, our DFT mo- 
ments are reduced ~10% from spin-orbit, but 50-100% 
by slightly reduced Fe-As spacing. Also, a Bogoliubov- 
de Gennes modeP^l suggests that only a large U induces 
domain walls. Yet, our results show that there are many 
low-energy magnetic defects, some with quenched mo- 
ments, even without invoking smaller Fe-As spacings. 

Nonetheless, spin fluctuations are agreed to play a big 
role, especially those involving magnetic planar defects, 

1. e., antiphase (APBs) and 90° domain (DBs) bound- 
aries (Figs. [I] and [2|. Such defects have been proposed 
as key to any moment description} 2 ^ but have yet to 
be quantitatively tested. With orthorhombic distortions 
(a > 6), structural and concomitant magnetic twins, see 
Fig. [21 are observed in BaFe2As2 along (110) with 10- 
50 \im between boundaries,^ or much less 100-400 nmW^ 
With stress samples detwin, but twins return upon its 
removal.^ As in YBa 2 Cu 3 07_^P^ (110) twins termi- 
nate on (110) twins. These defects cause anisotropic 



FIG. 1: (Color online) APB in the (a) be-, (b) ac- planes, and 
(c) 90° DB with no strain (a=b). Red (blue) circles are "up" 
("down") in-plane moments, as indicated in (a). HS, MS and 
LS indicate Fe-sites with high-, medium- and low-spin states. 



scattering near AFM wavevectors, giving 2-dimensional 
spin fluctuations. Twins also create stripes of increased 
diamagnetic response,^ and nucleated SC at magnetic 
defect boundaries.^ Niedziela et alW^ proposed that a 
high density of nanotwins in the a-b plane from local 
lattice distortions can improve the agreement between 
simulated and their measured pair-distribution functions 
(PDF). Local fluctuating, magnetic-structural variation 
introduces nanoscale AFM domains (coupling spin and 
lattice degrees of freedom) , potentially lowering the aver- 
age static moments. Thus, properties of magnetic planar 
defects (low-energy defects of the SDW), and related spa- 
tial/temporal averaging to estimate "ordered" moments, 
are critical to compare with scattering assessments. 

Here we use DFT to simulate various magnetic planar 
defects, i.e., two types of APBs, a 90° DB, twin bound- 
aries, and modified "nanotwin," which are low-energy ex- 
citations of the SDW. Figure and shows two APB 
boundaries in the Fe-plane - parallel to the be- or ac- 
planes - and an unstrained 90° DB is shown in Fig. [TJ^. 
Figure [2^i shows a typical example of an ideal twin. A 
modified nanotwin with 2-dimensional structural distor- 
tion (consistent with that suggested by Niedziela et alWty 
is shown in Fig. [2}d with a series of static displacements 
along a— and b— axis in the supercell. The undisplaced 
nanotwin with 1-layer of Fe separating defect planes is 
really a magnetic stacking fault (SF); a nanotwin super- 
cell has very different boundary conditions than a twin, 
with different far neighbors and distances between defect 
pairs; indeed, "ideal twin" supercells formed with 1-layer 
separation between defect planes (a high density of SFs) 
has local environments like the nanotwin, except that 
twin has symmetric relaxations governed by the super- 
cell periodicity, whereas the nanotwin has asymmetric, 
localized distortions to match the PDF. While we show 
the defect energies are similar, a nanotwin, due to its 
boundary condition and supercell, may be considered a 
fluctuating twin nuclei, which can have low-spin Fe-sites 
unavailable in the ideal twin supercell. 



FIG. 2: (Color online) (a) Twin boundary, without strain, 
separated by three Fe layers, (b) Nanotwin with boundary _L 
to (110), where atom positions in a distorted cell are barycen- 
tric weights of the cell corners. HS, MS and LS are indicated. 



For nondefected (parent) and defected cells we calcu- 
late energy per atom and the associated magnetic mo- 
ments (bulk moment is 1.6 /ie), and, hence, the planar 
defect energy, 7, defined as 7 = (i^def — Eo)d/V, where 
£?def and Eq are the total energy per atom of the defected 
and nondefected cell, respectively, d is the distance be- 
tween defect planes and V is the volume per atom. While 
the energy per atom is helpful, 7 is the appropriate com- 
parison for cost of creating the defect interface and its 
dependence on defect density and defect volume. 

To do this, we use VASP^with plane -wave pseudopo- 
tential projected augmented wave (PAW) basis J*^ with 
an energy cut-off of 380-420 eV. A Monkhorst-pack Bril- 
louin zone integration with a 16 3 k-mesh is used for the 
SDW (Fmmm) structure. Smaller k-meshes are used for 
supercells depending on the length along each axis. 

For APBs, we constructed doubled (2x1x1), quadru- 
pled, and octupled supercells to examine excitations, de- 
noted by 2- APB, 4- APB and 8- APB, respectively (Fig. 
[3). For APBs (Figs. [l}i,b), we use measured lattice 
parameters 1 (a=5.6146,_6=5.5742 and c=12.9453 A). For 
a 90° DB, we set a = b = (a + b)/2 = 5.5944 A to re- 
duce strain effects, and construct supercells similar to 
the APBs, denoted as 2-DB, 4-DB, and 8-DB. Twin 
(4[1 + n]x2xl) supercells (n=0, 1, ...) are denoted by 
(3 + 4n)-N Fe-layers between defect planes, and have 
4(1 + n) unit cells along a and 8(1 + n)xl0 atoms/cell. 
Nanotwin supercells are denoted 5-N, 9-N, and 13-N for 
Fe-layers between isolated nanotwin pairs; the supercells 
with the static displacements suggested by Niedziela et 
al. are more complex because the local distortions must 
be compensated within the cell (Fig. |2Jd). 

The energies and moments for APB and DB defects 
relative to the non-magnetic (NM) state are shown in 
Fig. [3] (top), and compared to the AFM ground state 
(SDW). In all cases, Fe moments have two behaviors: a 
high-spin state (HS in Fig. [T]) at sites away from bound- 
aries and a low-spin state (LS in Fig. [T]) at/near bound- 
aries. For APB (6c), the LS moment falls substantially 
to 0.8 jiB from 1.6 similar to that found by Yin et 



FIG. 3: (Color online) Energies relative to NM state (top) and 
Fe moments (bottom) for APBs (labels denned in text) in the 
(a) frc-plane (b) ac-plane, and (c) 90° DB. SDW indicates the 
AFM ground state. Insets depict the local Fe environment. 



FIG. 4: (Color online) Same as Fig. [3] but for (a) twins (la- 
bels defined in text), and nanotwins (3x3x2 supercell) with 
displacements along (b) da and (c) db axes. Fe has HS, MS, 
and LS states due to local environment (insets). 



a/P^ While for APB (ac), the LS moment decreases only 
to 1.54 jiB- Moments do not vary much with the size 
of the supercells, but these two structures energetically 
compete with the ground state SDW (< 9 raeV/atom). 
The two spin states depend on local magnetic environ- 
ments (inset Fig. [3}. For 90° DB (Fig. [J), the HS 
state has a higher moment of 1.7 /i# due to global strain 
from changed lattice parameters. The LS moment de- 
creases slightly to 1.57 \±b near the boundary. The local 
environment does not play a significant role, suggesting 
simple models such as counting the number of aligned 
neighbors is not sufficient to characterize the moments. 
This defect requires within 2 meV/atom excess energy to 
form compared to the SDW. They are energetically com- 
peting with the APB (ac). Both defects are then expected 
to be present at the same temperature. 

The energy and moments for twins are shown in Fig. 
|4^l. Interestingly, an Fe-atom in a twin has three spin 
states depending on the local environment. Fe- atoms at 
the boundary remain in a medium-spin state (MS in Fig. 
[2^i). A LS state occurs on Fe-sites adjacent to the bound- 
ary (Fig. [2|i). These Fe-sites have the same nearest- 
neighbor environment as the bulk HS states but differ in 
the farther neighbors. These defects can form at a few 
raeV/atom, albeit 7 is more critical, see below. 

The nanotwin energies and moments versus distortion 
along a— and b— axis (in A unit) are shown in Fig. j4^b,c). 
Similar to twins, there are three Fe spin states: a HS bulk 
state (1.6 /ig), a MS state (1.42 /i^) at the boundary, and 
a LS state (0.8 — 1.0 /i#) in the vicinity of the distorted 
side of the boundary. The structural perturbations show 
a stronger effect on the LS moments near defect bound- 
aries, decreasing to as low as 0.8 Isolated (fluctuat- 
ing) nanotwins are equally competitive to form as dense 
twins but with much reduced moments. Energies are 
affected mostly by the changed magnetic configurations 



and very little by distortions. So, magnetic defects drive 
the short-range structural distortion (not the other way 
around) and can help quench magnetization. 

Table [I] compares the various planar defect energies. 
Small APB-ac, 90°-DB, and isolated nanotwin defects 
(low-energy spin excitations) compete with widely sepa- 
rated twins (spin kinks), which are observed. Energeti- 
cally, APB (ac) and 90° DB are the most favorable and 
most likely to persist after annealing, but twins form and 
are stabilized by lattice strain from disclinations from 
twins oriented in the two independent directions. 33 Such 
small fluctuating defects will affect the observed average 
moments, whereas separated twins will affect the mag- 
netic correlation length. Planar defect energies (7 in Ta- 
ble give the relative order in which magnetic defects 
can form and remain after processing. Structural defects 
can act as pinning sites for magnetic domain walls. 

Strikingly, separated twins are higher in energy than 
dense twins, until a separation of 16 unit cells (15 Fe- 
layers), where 7t w in reaches a maximum (Table [i]), after 
which there is a slow convergence of 7t w in versus d (Fig. [5] 
and 00]). At 28 unit cells (~11 nra), 7twin has not yet 
converged; at best convergence could be reached at 32 
unit cells (~13 nm), more likely larger - depending on 
if 7 drops near that of densely packed twins. Observed 
structural twin^MiLU are extended well beyond the ones 
computationally possible. Twin separation L oV is affected 
by stress, where equilibrium is reached with applied stress 
roughly 7t w in^- Also, in real samples twins appear in 90° 
oriented pairs, reducing strain at the cost of increased L d\ 
While separation will depend on local defects and stress, 
the average distance L oV between magnetic twins (esti- 
mated from Fig. [5] at 32+ unit cells) should yield a peak 
in the magnetic susceptibility. Indeed, the measured 
magnetic correlation length is ~30 unit cells ! 3Q * 31 IWe sug- 
gest that small, low-energy excitation can reduce average 



4 



TABLE I: 7 (27 for twins) for various planar defects (in 
mJ/m 2 ). Energies (raeV/atom) are relative to SDW. 7twin 
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FIG. 5: Twin energy versus number of Fe- layers separating 
twins, or distance between twin planes. 

moments by spatial and temporal averaging, while the 
twins dictate the magnetic correlation length. 

A nanotwin (Fig. ^p) with no local structural dis- 
tortion is like an isolated, ideal defect pair, not a dense 
set of twins. To understand the effect of short-ranged 
structural distortion, we have studied a 1-N ideal twin 



with(out) relaxation in afr-plane for only those atoms 
near the boundary, more localized than in the nanotwin 
supercell. The planar defect energy with(out) relaxation 
is 62 (74) mJ/m 2 . The relaxation along a- and 6-axis 
lie within 0.9% of the ideal case, close to the best fit 
to measured PDF, 36 so the twin and nanotwin are very 
similar in energy and local structure. However, a nan- 
otwin may be considered a fluctuating twin nuclei, which 
has many more LS sites (Fig. [4| not available in a twin 
supercell. Unlike for ideal twins, the nanotwin surface 
energy decreases to its limiting value as the nanotwin- 
nanotwin distance grows. For insight into an electronic 
origin of these magnetic effects we have calculated the 
Fe-site projected density of states (see [40]). 

In summary, we have studied competing low-energy 
magnetic planar defects in BaFe 2 As 2 . The favorable de- 
fects are APB(ac), 90°DB, and nanotwins, but twins (as 
observed) are favorable on a longer length scale. The 
most pronounced reductions in Fe-moment are near the 
boundaries of APBs (be) and dense twins. We find that 
isolated closely-spaced twins are energetically favorable 
and correspond to a recently proposed 36 nanotwin sug- 
gested to match the PDF from experiment. The low- 
energy defects in the perfectly ordered SDW, and the as- 
sociated spatial and temporal averaging to estimate the 
"ordered" moments, are critical for comparison to exper- 
imental estimates, and serve as a possible explanation 
for the discrepancy between the Fe moments predicted 
by DFT (1.6 - 1.9 fi B ) and experiment (0.8 - 2.1 fi B ). 

In DFT there is no global quenching of Fe- moments, 
rather quenching is localized near defect boundaries, sup- 
porting the increased diamagnetic response^ at bound- 
aries. Four things can reconcile the various results: (i) 
the sample may have a high density of magnetic nan- 
otwins that yield strong reduction of the average mo- 
ment; (ii) the low-energy magnetic defects may be weakly 
mobile, affecting any time averaged value; (iii) magnet 
defects with fluctuations in the local Fe-As planar spac- 
ing impact values of local moment; and/or (iv) the low- 
energy fluctuations around the SDW have an extended 
energy window, which is not currently used in scatter- 
ing experiments to estimate local moments from inte- 
grating the susceptibilities. Indeed, the local moment 
in CaFe 2 As 2 (0.5 fi B ) 41 and BaFe 2 As 2 (0.5-0.8 /i B )P 
were estimated by integrating up to about 0.1 eV, where 
DMFT suggests ~1 eV window. 42 It would be beneficial 
to adjust the energy window in neutron-scattering to as- 
sess this effect on the assessed Fe moments. 
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